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We discuss possibilities of measurement of deeply virtual Compton scattering ampli- 
tudes via different asymmetries in order to access the underlying skewed parton distribu- 
tions. Perturbative one-loop coefficient functions and two-loop evolution kernels, calcu- 
f^ • lated recently by a tentative use of residual conformal symmetry of QCD, are used for a 

^sj I model dependent numerical estimation of scattering amplitudes. 

1,^ . 1 Introduction. 

vQ I The asymptotic freedom of QCD allows to calculate short distance quantities 

^+ ■ perturbatively form the first principles. Together with the use of factorization the- 

^^ I orems, the cross sections for (semi-) inclusive reactions are predicted at a large scale 

as a convolution of perturbative parton cross sections and long distance parton 
distributions. In the case of exclusive reactions, such as two-photon reactions in 
the light-cone dominated region or diffractive processes, the amplitudes factorizes 
i-C I in hard-scattering amplitudes and hadron distribution amplitudes and/or skewed 

parton distributions (SPD). Such distributions were already proposed in the past 
Q^l [Q and they become more popular nowadays [g] since it has been realized their 

1^ ■ vital role for the proper description of the mentioned reactions. Moreover, as it was 

pointed out by J. Xi, their measurement can clarify the issue of how the hadron 
total angular momentum is distributed among its constiutents. One can expect that 
the study of the SPDs will be one of the central issues of hadronic physics in nearest 
future. 
C^ I A promising way to access the SPDs is to probe the nuclcon with a highly virtual 

photon via the Compton scattering process, which is called deeply virtual Compton 
scattering (DVCS). Indeed, it is auspicious that a first evidence of experimental 
feasibility to detect such a signal has been revealed by the ZEUS collaboration 
Pi . On the theoretical side the factorization of the coUinear singularities has been 
proved at leading twist- two level to all orders in coupling El. The local version of 
the OPE does not converge in the considered kinematical region. Fortunately, this 
problem can be overcome by a non-local light-cone expansion in terms of light- 
ray operators 0, W. Moreover, from the renormalization group equation of these 
operators one can derive in a straightforward manner the evolution equations for 
SPDs. 
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Here we shortly review our approach to complete next-to-leading (NLO) cal- 
culation, which is accomplished by a heavy use of the conformal symmetry to 
map the known results from forward deep inelastic scattering (DIS) to the non- 
forward kinematics of DVCS and simple one-loop addenda. The problems within 
our modus operandi have been resolved by employing conformal Ward identities 
and constraints arising from the algebra of the coUinear conformal group H . 



2 The virtual Compton scattering process. 

The virtual Compton scattering (VCS) process on a nucleon can be measured 
in the reaction 

£^N -4 £^Nj. (1) 

However, as shown in Fig. y, there is contamination from the so-called Bethe-Heitler 
(BH) process in which the final real photon is radiated from the lepton beam. Thus, 
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Fig. 1. The virtual Compton scattering process (a) and the Bethe-Heitler process (b). 
the scattering amplitude square jTp in the cross section. 



contains beside the VCS and BH parts also the interference term: 

\T\ — 2_^ \\Tvcs\ + |7s_ff I ^TvcsTbh^T^vcs^bh) 



(3) 



The BH-amplitude is purely real and is given as a contraction of the leptonic tensor 

L^, - S(fc', A') [7^(1^ -l^Yhv + 7.(r +f )"Sp] "(fc, A) (4) 

with the hadronic current J^ = u^P^) {Fx{^^)^v + iF-2,{^^)ayr^^ u{Pi): 



TBH = --r^e*(A)^L^'^J., where A = Pi - P2 



(5) 
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Since Dirac and Pauli form factors Fi and F2 are (partially) known from exper- 
imental measurements and can be parametrized to a certain accuracy by double 
dipole formulae, one can disentangle the content of the hadronic tensor Tycs'- 



f? .... ^..u^.nj^ '- ^ f- for 



rycs = T-e*(A)^T^''u(fc')7.w(fc), where <^ , ^ , . (6) 

qi |_ + tor e+ 

This tensor, defined by the time ordered product of two electromagnetic currents 



contains for a spin-1/2 target twelve^ independent kinematical structures H. 

Measuring the cross section (El) in quite different experiments: with negatively or 
positively charged lepton beam, with (un)polarized beam or (un)polarized target, 
the rich structure of the hadronic tensor can be explored in different kinematical 
regions. In the following we are interested in DVCS, where —qf is large and A^ 
is small. Unfortunately, since the DVCS amplitude (O) decreases with increasing 
qf while the BH amplitude (0) grows up with decreasing A^, in this kinematics it 
will be a delicate task to extract the VCS amplitude for large virtualities —qf and 
especially low A^ from experimental data. On the other hand it is promising to 
look for asymmetries which allow one to measure the interference term directly and 
get in this way access to the real and imaginary part of the DVCS amplitude [pi: 

1. charge asymmetry 

da{e~-hea.m) — dcr(e^-beam), oc 4ReTDvcsTBH (8) 

2. single spin asymmetries (either the lepton-beam or the nucleon is polarized) 

da{t) - da{i) ex 4ImrDvcsrBH, (9) 

3. azimuthal angle asymmetry 

f ,„ d(hrda — f ,„ dSrda 

A = ^-^1^ ,^ -^"/' , (10) 

J„ d(prda 

where (f)r — 4'N~4'e with (f)]y{(l)e) is the azimuthal angle of the final nucleon (lepton). 



^) 12 = i X 3 (virtual photon) X 2 (final photon) X 2 (initial nucleon) X 2 (final nucleon). The 
reduction factor 1/2 is a result of parity invariance. 
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3 Practical implication of conformal symmetry. 

To describe different "two-photon" processes in the light-cone dominated region 
in a unique way, we generalize the Bjorken limit and introduce the scaling variables: 

^^1^, ry^^, Q2^-g2^-i(,i+,,)2, P = P,+P,^ (11) 

where ^ « X'q-J{2 — xbj) is related to the Bjorken variable xbj — —q\/2{Piqi) 
and f] is the so-called skewedness parameter. Formally, this parameter interpolates 
between different processes as given in table 0. 



DIS 
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V = 


= 




DVCS 




q',=0 


V = 


= l/uj 
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Table 1. The value of the scaling variable rj for different two-photon processes. 

Since in the case of DVCS the longitudinal part of the hadronic tensor vanishes, 
there are only two kinematical structures left to leading twist-two accuracy: 

T^u = -gl,:Fi{^,v = ^,Q^A') + ^e^.gp^Qi{^,v = ^,0", A^) + • ... (12) 

Here the transverse part of the metric tensor, denoted by gj^^, and the e-tensor 
are contracted with the projection operators "Pq^ = ga/s — qiaq2i3 / {11(12) '■ 9^v — 
'^/i"5Q/3^1' ^tJ.i^qP = 'PfT^a/Bqp'P^,,- The appearing amphtudes can be calculated by 
means of the OPE. In contrast to the case of DIS, a larger number of (local) oper- 
ators and their Wilson coefficients contribute in the non-forward case. Fortunately, 
we can appreciate the conformal OPE (COPE) Q, which tells us that all Wilson 
coefficients are fixed up to an overall normalization, which can be taken from known 
results in DIS. The conformal operators appearing in the COPE are: 



QqV 

QqA 



^i^9+y\j+_ ^C/I^U, (13) 



GqV^ 

gqA> = G+f,{td+) 




where d =d + d and D=D — D- The + and — components are obtained by 
contraction with the two light-like vectors n and n*, such that n^ = n*^ = and 
nn* = 1. They carry the spin quantum number I + 1 and the conformal spin j + I. 
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In the conformal subtraction (CS) scheme their mixing is only induced by the trace 
anomaly and proportional to (3 llQ] : 

/^^<^.'=7,,C?,, + ^E^.''^^fc' (14) 

The anomalous dimension matrix 7 is, up to a normalization factor in the mixed 
channels, the same as in DIS. 

Unfortunately, for the case of DVCS kinematics the COPE in its local version 
does not converge. Thus, it is necessary to resum the local operators into so-called 
light-ray ones. This provide us a definition for the SPDs: 

X|(i,^) = /^e^'^*^+(P252|^(-«n)| ^+ \tP{^n)\P,Si), (15) 

^q J J 27r 1^ 7+75 J 



t+V p tt-ri 



G„Af (*''?) = — y'ge-*^-(P252|G%(-Kn)|^^JJ^"^^|Ge+M|Pi5i). 



The momentum fractions of the incoming [outgoing] parton is j^P^ [j-tt-P+J. A 

form factor decomposition would give us the functions H, E and H , E for the 
parity even and odd sectors, respectively, as introduced by Ji [Q. 
Next we have to reconstruct the coefficient functions appearing in 



Q 



NfV 



(16) 



where T^ — J-'i, T — Qi and the sum runs over quark species Q — u,d,s with 
electrical charge eg. The convolution is defined by ri ® T2 = ti{x, y)T2{y, z). This 
task can be achieved by some manipulation of the COPE and, e.g. for the flavour 
non-singlet case, can be cast into the form 



r(f,r;,t,Q^M') = 7^ -'- ®pT ®C -,-;aJ. (17) 



The evolution kernel V (-, -) and the coefficient function C (f, -;as] are diagonal 
w.r.t. the conformal waves — the Gegenbauer polynomials: 

-^'^e)«{c}a.::-^{-:i:!H^i::). <'«> 

and their eigenvalues 7j(as) and Cj{as) coincide with anomalous dimensions and 
Wilson coefficients, respectively, which appear in DIS. A closed form of the function 
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F can be deduced order by order in as from the sum 



F(x,2/) 
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(19) 
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Conformal covariance does not hold true in the minimal subtraction (MS) 
scheme. The transformation, which rotates the CS result (|l^ -|l^) to the mini- 
mal one, is induced by the conformal anomalies appearing in the MS scheme. The 
rotation matrix can be found by study of the conformal Ward identities for Green 
functions with operator insertions and calculating the anomalous contributions via 
(modified) Feynman rules. Indeed, the off-diagonal piece of the anomalous dimen- 
sion matrix to NLO is induced by the special conformal anomaly matrix 'yC^''^ and 
the first coefficient of the /3- function m : 



1 



(1) - -^,(l) 
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, b,k = - [1 + (-1)-'"1 ^,fc(2fc + 3)1. (20) 



The remaining task to transform the coefficient functions to the MS scheme 
have been solved already in ||l^. In the case of DVCS, they are nothing else as the 
analytical continuation of the hard-scattering amplitudes for the meson production 
by two photon-collision. In the case of the axial-vector channel they read: 
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(21) 
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4 Numerical estimations. 

There are different proposals for modeling the SPD. Here we will deal only 
with two of them, equating the SPD equal to the forward parton distributions 
(FPD-model) and using the proposal of Radyushkin for the ansatz of the double 
distribution (DD) functions ||ll| /(z_,z+) = 7r(z_, z+)/(z-)_), where the profile 
function tt for quarks and gluons are given by 



3 3 [l-|z+|]^- 

"("-'"+) = 4 [l-|z+| 



7r(z_,z+) = — ■ -T , (^11) 



16 



[i-k+|] 
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The SPD is obtained by an integration M : 

.1 i-l-\z+\ 

q{t,v,Q^)= dz+ dz^S{z+ + r,z^-t)fiz^,z+,Q^). (23) 

J-l J-l+\z+\ 

We apply the DD model at different scales. To evolve the SPD^ we expand them: 

(24) 



3=0 



in terms of Legendre polynomials and express the expansion coefficients Cj{r],Q'^) 
in terms of conformal moments which evolution is governed by the anomalous 
dimension matrix (EQ). This method is suited for not too small 77. 







Mim 0.75 




Fig. 2. In (a) and (b) we show the SPD for the sum of the non-polarized u and u quarks 
and the gluon densities, respectively. The FPD-model (dotted line) and the DD-model 
with MRS llSl parametrization (dashed line) are modeled at a scale of 4 GeV , while the 
DD-Model with GRV H] parametrization (solid line) has been taken at the momentum 
scale 0.4 GeV^ and has been evolved afterwards with NLO formulae to Qq = 4 GeV'^. 
The skewedness parameter is ?7 = 1/3 which corresponds to xbj = 1/2. 



qf = AGeV^. 



In Fig. we show our input distributions at a scale of Q^ = 
Let us also mention, that the relative NLO corrections due to the evolution are as 
expected not very large. Starting at an input scale of 0.7 GeV, we found that the 
NLO corrections do not exceed 10% to 30% ||. 

For the numerical estimation of the hard-scattering amplitudes ( |l6| ) we perform 
also an expansion with respect to Legendre polynomials which allows us to include 
easily evolution effects. However, this method is not appropriate to study the low 
XBj region. In general we found that the NLO corrections for a;Bj > 0.1 are typically 
of the order of 20-50% (see Fig. ^). In special cases, e.g. for the DD model with 
MRS parameterization, it can even be much larger. Fortunately, there is partial 
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Fig. 3. The real part (a) of T\ for DD-model with MRS and GRV parametrizations for 
the scale Qo = 4 GeV^ and for xb.\ > 0.1. Here MRS (GRV) amplitude is plotted at LO 
as dashed (dotted) and at NLO as solid (dash-dotted) lines, respectively. The ratio of the 
real to imaginary part (a) of T\ , respectively, for the DD- and FPD-model with the MRS 
parametrization. Here the DD(FPD)-model at LO and NLO is shown as dashed (dotted) 
and solid (dash-dotted) lines, respectively. For the SPD-models considered the result is 
independent of A'^. 



cancellation of perturbative corrections in the ratio of real to imaginary part. As 
shown in Fig. Hb, this ratio is model dependent and may its measurement can be 
used to get a deeper insight in the SPD. 
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